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Abstract 



We analyze the structure of the periodic trajectories of the matrix gen- 
erator of pseudorandom numbers which has been earlier proposed in 0,0]. 
The structure of the periodic trajectories becomes more transparent when the 
rational sublattice coincides with the Galois field GF[p] || ||]. We are able 
to compute the period of the trajectories as a function of p and the dimension 
of the matrix d. 
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1. In the articles |2| the authors suggested the matrix generator of pseu- 
dorandom numbers based on Kolmogorov-Anosov K systems ||, f| [7[ ||. These 
systems are the most stochastic dynamical systems, with nonzero Kolmogorov en- 
tropy H |, [11], [L2], The trajectories of the system are exponentially unstable 
and uniformly fill the phase space. The coordinates of these trajectories represent 
the desired sequence of pseudorandom numbers fll], 0. 

In the given case it is assumed that the unit d- dimensional torus Tl d plays the 
role of the phase space, therefore the coordinates of the trajectory are uniformly 
distributed over the unit hypercube and can be used for Monte-Karo simulations 

[0, 0]- 

The properties of this new class of matrix generators were investigated with 
different criterion including discrepancy in various dimensions. In all cases it 
shows good statistical properties M. Matrix generators based on different ideas are 



proposed in (TTJ], [TJJ. We refer to the book of Niederreiter []TB| and to the survey 
article |L7| for recent references. 

The aim of this article is to estimate the period of the trajectories which are 
used to produce pseudorandom numbers. Only periodic trajectories of the system 
can be simulated on a computer, because trajectories on a computer are always on 
a finite rational sublattice Z p of the phase space U d |L0], [LJ 



The essence of the approach is to consider the system on rational sublattice 
Z p of a unit d-dimensional torus and particularly on sublattices with prime basis p 



20fl . These sublattices are equivalent to Galois fields GF[p] and all four elementary 



arithmetical operations can be carried out unrestrictevely [F2T], |23|]. This approach 
demonstrates that in order to have a large period on a sublattice GF[p] one should 
have matrices A which have the eigenvalues in high extensions GF[j¥p] of the field 
and at the same time they should have integer entries. 

On every Galois field exists a primitive element g which has the period p — 
1, that is g v ~ x = 1, and every element of GF\p] can be represented as a power of 
g [[H], |2^|. We shall show that if the eigenvalues of the matrix generator are 
proportional to these primitive elements, then the period of the trajectories is equal 
to 



r p (A 2 )=p-l. 

The upper bound for the period of the cat map (15) found by Dyson and Falk 
on arbitrary rational sublattice Z q is equal to 



r <3q 

and also linearly depends on q. 

The period increases when we consider the sublattices Z pn of the basis p n , where n > 
2 instead of p. In that case the period of the same matrix generator is 

T pn (A 2 ) = (p-l)p n ~\ 

Then increasing the dimension of the matrix generator by {2d x 2d} matrices we 
will get the period equal to 



r p {A 2d ) =d(p-l) 



and correspondingly Oil Zj<nn sublattice 



T P *(A 2d ) = d(p-l)p n -\ 

Using quadratic extension of the Galois field GF[^fp\ |21|, ^2], |23| we have found a 
systematic way to construct the matrix generators with period 

t p = p 2 — 1, and T p n = (p 2 — l)p n ~ l 
and finally on the high extended Galois field GF[^fp] the period is equal to 

t p = p d - 1, and Tpn = (p d - l)p™ -1 . 

The last result shows that in practical simulations the period is very large and is of 
order p d+n ~ l where d- is the dimension of the matrix generator andp n - is the basis 
of the sublattice Z p n . 

We suggest specific matrices with this properties and with almost zero entries, 
see Section 17. This matrices have the largest period and can be easily used for 
practical simulations. 

2. Let us pass to the details of the algorithm. The matrix generator is defined 
as§|, 

X (n+1) = A ■ X {n \ (modi), (1) 

where A is d x d dimensional matrix with integer matrix elements a^- and determi- 
nant equal to one 

Det A — 1, (2) 

and = (X{ 0) ,...,xf ) is an initial vector. The last condition provides phase 
volume conservation. The automorphism (1) form the K system of Anosov if and 
only if all eigenvalues of the matrix A are in modulus different from unity || |7j, [§| 

I Ai| 7^1, i = l, -,d (3) 
The trajectory of the K system (1) 

represents the desired sequence of the pseudorandom numbers [JE]. 

This approach allows a large freedom in choosing the matrices A for the K system 
generators and the initial vectors [0. Specific choices suggested in |], @, ^ are 



( 2,3,4, ,d ,1 \ 

1,2,3, ,d- 1,1 

1,1,2, ,d-2,l 



1,1,1, ..,2, 3, 4,1 
1,1,1, ..,1,2,2,1 
1,1,1, ..,1,1,2,1 
V 1,1,1,.., 1,1,1,1 / 



,A d = 



( 0, 1, 0, , 

o, o , 1 , , 



0, , , , 1 

V (-l) d+1 ,ai,a 2 ,.,a d _i J 



(4) 



3. Let us consider trajectories of the system (1) with initial vector which 
has rational coordinates 

Xo = ( ^ K..,^ )• (5) 

Pi P2 Pd 

It is easy to see, that all these trajectories are periodic orbits of the Anosov map 
(1), because matrix elements ajj are integer. Indeed, if we shall consider the 
sublattice of unit torus U d with rational coordinates of the form q/p where p is 
the smallest common denominator 

X = ( ^ q A-^ ), 0< qt <p-l (6) 
p p p 

then the multiplication, summation and (mod) operations (1) will leave the trajectory 
on the same sublattice. The total number of vertices on this sublattice Z d is 

(total number of verteces) = p d , (7) 

therefore the period t p of the trajectories on Z p = Z p <g> ... <g> Z p , where Z p = 
{0, 1, ...,p — 1} is always less than p d 

r p < p d . 

Thus the periodic trajectories of this system (1) with the initial vector (5) coincide 
with a subset of the points of rational sublattice Z d and our goal is to find conditions 
under which the period of the system will be as large as possible. 

Let us show that on every given sublattice Z d Anosov map (1) reduces to (mod 
p) arithmetic. Indeed on sublattice Z d the Anosov map A (1) can be written as 

(n+l) in) 

^ = E«m (mod 1) (8) 

and is equivalent to (mod p) arithmetic on the lattice with integer coordinates qi 
which are in the interval [0,p — 1] 

# +1) = E«M?! n) , (mod p). (9) 



Thus the images of the periodic trajectories on a unit torus U d appear as trajectories 
on the integer sublattice Z d and all operations can be understood (mod p). The 
most important thing is that now all operations become commutative. 

4. To estimate the period of the trajectories on rational sublattice it is essential 
to consider those sublattices for which p is the prime number, we mean that p\ = 
... = p d = p. In that case the integer sublattice gains an additional structure and 
becomes the Galois field and all operations reduce to arithmetic ones on Galois field. 
The benefit to work on Galois field is that four arithmetic operations are well defined 
on that sublattice [21] . 



In this way we can consider every coordinate , i — 1, ...,d as belonging to 
Galois field GF[p] = {0, 1, p — 1}, where p is prime number and consider the 
sublattice direct product of Galois fields 

Z$ = GF\p) ® ... ® GF\p]. (10) 

As we already mentioned, this reduction of the dynamical system (1) to dynam- 
ical system for which the Galois field plays the role of the phase space makes all 
operations commutative in the sense that 

{A{A X}} = {A 2 X}, (11) 

where {...} means mod operation. The commutativity of the multiplication and 
(mod) operation on the Galois sublattice means that the periodic trajectory 

WA {AX}...}} = X (12) 

can be represented in the form 

{A Tp X} = X (13) 

and the period of the trajectory r p can be understood as a degree of power on which 
the matrix A reduces to identity (mod p) 

A Tp = 1 (mod p) (14) 

The period of the trajectory on the Galois sublattice is equal therefore to the power 
r p in which the matrix A reduces to identity on a given Galois field GF[p\. This 
period does not depend on initial vectors and the whole phase space Z d factories 
into trajectories with equal periods. It is obvious that the same matrix A will have 
different periods on different Galois fields and that this period depends on the given 
prime number p and the dimension d of matrices. 

5. To demonstrate this fact let us consider few examples. The matrix 

^=(?'JV ( 15 ) 



has period equal to four on the Galois field with p = 3 and to eight when p = 7 



r 3 (A)=4, r 7 (A) = 8. 



The question which appears here, is how it is possible to estimate the period of the 
matrix A without actual computation of the powers of the matrix A. 

We can find the answer to this question considering the eigenvalues of the matrix 
A. Indeed, as we will see, we can compute the periods using the eigenvalues of the 
matrix A. Let us consider first the example of the eigenvalues of the cat map (15) 

. 3 + 3-^5 
+ 2 ' 2 
The question: what is the period of the given matrix A on Galois field GF[p] is 
equivalent now to the question : in which power the eigenvalues are equal to identity 
on field GF\p\? 

\ t p = 1 (mod p) (16) 

As it is easy to see 

A+ = 1 (mod 3), A+ = 1 (mod 7) 

which confirms the direct computation. The exceptional case when this method can 
not be applied directly is when the eigenvalues have degeneracy on a particular field 
GF[p\. This takes place for p — 5, indeed A + = A_ = 4, (mod 5). Using 
Jordan normal form of the matrix one can see that r 5 = 10. Because this happens 
for very particular values of p in the following we will consider only the cases when 
eigenvalues are not degenerate. 

6. Thus the actual value of the period r p naturally depends on the form of 
eigenvalues A and of the prime number p. Here we can distinguish different cases: 

i) . The eigenvalue A coincides with one of the elements of the Galois field GF[p\. 
In that case the period r p depends on whether eigenvalue coincides with the primitive 
element of the Galois field or not. All elements of the field GF\p\ can be constructed 
as powers of primitive element g and g p_1 = 1. If one of the eigenvalues coincides 
with the primitive element of the Galois field , 

A = g, where g is the primitive element of GF\p], (17) 
then the period of the matrix is maximal and is equal to r p = p — 1 

A p ~ x = 1, (mod p). (18) 

Therefore to get the maximal period in the case i) one should have at least one 
of the eigenvalues equal to the primitive element of the field GF[p\. If A does not 
coincide with the primitive element g, then the period is simply smaller. 

ii) . The eigenvalue does not coincide with any of the elements of the Galois field 
GF[p\. This takes place because the solutions of the characteristic polynomial of the 



matrix A are not always in the field GF[p\. Galois field is arithmetically complete, 
but it is not algebraically complete, therefore one can have the situation when 

A is not an element of GF[p\. (19) 

This possibility can be illustrated by cat map (15), indeed a/5 is not an element of 
GF[3] or GF[7]. 

In that case one should ask, whether it is an element of the quadratic extension 
GF[y/p\. The quadratic extension of the Galois field consists of the numbers of the 
form a + b^fg where a, b are the elements of field GF\p], g is the primitive element 
of GF[p] and yj~g is a square-free integer. 

Now if the eigenvalue is an element of the quadratic extension and coincides with 
it's primitive element h 

A = h, where h = hi + h 2 \fg is the primitive element of GF[- s /p\, (20) 
then the period is equal to r p = p 2 — 1 

X? 2 - 1 = 1, {mod p) (21) 

because the primitive element of the GF[y/p] has the period equal to p 2 — 1 pl| . 
Again, if A does not coincide with the primitive element h, then the period is simply 
smaller, as it is in the case (15) for p = 3, 5 where the period is p + 1. 

iii). In general the characteristic polynomial of the matrix A is of order d and the 
eigenvalues can belong to high extensions GF[^/p] of the Galois field, the elements of 
which have the form a + bq + . . . + eq d ~ x where a, b, . . . , e are the elements of GF [p] , g 
is the primitive element and q d = g. The primitive element h = hi + h 2 q + ... + h c iq d ~ 1 
of GF[j¥p] has the period r p = p d — 1 

A^ 1 = 1, (mod p). (22) 



This analysis demonstrates an important fact that in order to have a large period 
on a sublattice GF[p] one should have matrices A which have the eigenvalues in 
high extensions GF[tfp] of the field and at the same time they should have integer 
entries. 

7. In the previous sections we described the trajectories of the K system on the 
rational sublattice Z d and particularly on a Galois field, that is when p is the prime 
number. 

In this section we will reverse the discussion and will try to construct the matrices 
A with the properties of K systems on a given Galois field with the maximal period. 
The question can be formulated in the following form: can we construct a matrix A 
with the properties of K system such that it has the largest period on a given Galois 
field GF[p] ? 



Let us first consider two-dimensional matrices of the form 

^(Virt 1 )- (modp) (23) 

which have the following eigenvalues 

A + = a + Vex 2 — 1, A_ = a — Vet 2 — 1 

To realize the first case i), when the eigenvalue belongs to the field, we should have 
Va 2 — 1 as an element of the field, that is 

a 2 - 1 = k 2 , k^O. (24) 

In this case the square root operation will belongs to the field. To have the maximal 
period we should choose one of the eigenvalues to be the primitive element g of the 
given field GF[p] 

X + =g, X.=g-, (25) 

therefore 



a + k = g , a — k = g 
and (24) is satisfied automatically. From (26) 



(26) 



9 ~9 



9+9 , - 

« = — ^ — , 9 T 9 



(27) 



so that the matrix (23) is equal to 



.4, 



/ 9~+9 
2 



9 ±R + 1 \ 

, 2 



V 



9 +9 _ 1 9 +9 

2 ' 2 



(28) 



/ 



and has the period as large as the primitive element g, which is p — 1 

r p {A 2 )=p-l. (29) 

The upper bound for the period of the cat map (15) found by Dyson and Falk ]nj 
on an arbitrary rational sublattice Z q is equal to 

t <3q 

and also linearly depends on q. 



8. These formulas allow to construct explicit examples of matrices with given 
period. The field GF[7] has primitive element g = 3 and g" = 5 therefore the matrix 
(28) has the form 



A = l Z ° A , (mod 7) (30) 




with the corresponding period r 7 = 6. 



9. The next step in this construction is to enlarge the sublattice GF[p] to sub- 
lattice Z p n, where p is the same prime number. Despite the fact that the sublattice 
Z p n does not have field structure, nevertheless there exists an element h with the 
period ( p — l)p n ~ l . The important theorem |21[] states that h coincides with one 
of the primitive elements of the original field GF[p] which has the property 

ST 1 ^ 1, {mod p 2 ). (31) 
This primitive element g is the same for Z p « for any integer n > 2 and has the period 

We have therefore the following result: the matrices which we have constructed in 
the previous section on GF[p], will have period on sublattice Z pn equal to 

r pn (A 2 ) = (p-l)p n - 1 (mod p n ). (32) 

It simply means that with (mod p n ) operation we increase the period of the matrix 
A 2 from p — 1 to (p — l)p n_1 . This allows to have large sublattices with small basic 
primes. 



10. For the case GF[7] the condition (31) is satisfied because 

3V 1, (mod 49) 
and the matrix (30) has the period 

r 7n (A 2 ) = 6 • T- 1 for anyn>2. 

So our construction of the matrices with eigenvalues which are proportional to the 
primitive element of the field GF[p] is completed. 



11. As a basis for next constructions, let us consider a class of K system 



generators with very simple structure 24 



/ 0, 1, 0, 
0, , 1 , 



A,, 



0, , , 
w+i 



o \ 





(33) 



V (-l) d+1 ,a 1 ,a 2 , ..,a d -i J 
In the last case it is easy to compute the characteristic polynomial of Ad 



A d — a<i-i A 



d-i 



ai A + (-l) d = 



(34) 



and therefore for it's eigenvalues Ai, Xd we have 

Ai • • • X d = 1 



Ai + ... + A d = a d _i. (35) 

These formulas allow to choose eigenvalues and then to construct matrix Ad for K 
system generators. 

For example if d=4 and a\ = 0, a 2 = 3, and = 0, then 



2 

with an additional simplectic structure of A 4 . 



As = , i^I, A, = - Ji^, (36 ) 



12. Our goal is to get matrices of the form (33) with the maximal period. Let 
us first consider four-dimensional case 



( 0, 1 , 
0, , 1 
0, , 





1 



(37) 



V -1, 0, 2a, / 
which has the characteristic polynomial 

A 4 - 2a\ 2 + 1 = 0, 
and we will choose again 2a = g + g~ , then the roots are: 



(38) 



Ai, 2 = ±y/g, A 3 ,4 = ±V<T- (39) 
Because g is square-free primitive element of the GF[p] the period of this matrix is 



r p = 2(p-l). 

Increasing the dimension of the matrix A with the same simplectic structure 

\ 



(40) 



I o, 1, 0, 
0, 0, 1, 



i2d 








(41) 



0, , , , 1 

V -1, , ,...2a,.., J 

we will get the characteristic polynomial 



\ 2d - 2a\ d + 1 = 0, 2a = g + g- 



(42) 



with 2d different roots 

X d + = a + Va 2 - 1, \t = a- Va 2 - 1 (43) 
and the period is equal to 

r p (A 2d )=d(p-l). (44) 
The same matrices A 2 d on Z p n sublattice will have the period 

T pn {A 2d )=d{p-l)p n -\ (45) 



13. The example on GF[7] where g + g~ 



( 



A 2d = 



0, 1, 0, 
0, 0, 1, 








8 



1 = — 6 is 



(mod 7) 



o, o , , , 1 

V -1, , ,..-6,.., J 

with period r 7 = d ■ 6 and with (mod 7 n ) we have r 7 n = • 6 • 7 n_1 . 



(46) 



14. The next step is to construct the matrices which have the eigenvalues in 
quadratic extension GF[y/p], that is we are going to consider the case ii). If h is the 
primitive element of the GF[y/p], that is 

h = h 1 + h 2y /g, h-h* = g, h + h* = 2h u (47) 

then the matrix which has the eigenvalues in GF[^/p] can be constructed in the 
same form as (33) 



0, 1, 

0, 0, 1 

\ -1, 2h ig --g, 2h x -g- J 

because the characteristic equation is 



(mod p) 



(48) 



A 



(\ + g-)(\-h)(\-h*) = 

(2hi - g-)\ 2 - (2h 1 g~ - g)\ + 1 = (mod p) (49) 

and has the root h which coincides with the primitive element of GF[y/p\. This 
matrix has integer elements by construction and the period 

r p (A 3 )=p 2 -l. (50) 

Therefore the period quadratically increases in comparison with previous construc- 
tion. The same matrix with (mod p n ) operation will give 



r pn (A 3 ) = (p 2 - l)p n ~ 



(51) 



15. The example on GF[7] where h = 1 + 2 ■ v3 will be 

(mod 7) 



A, 



1- 
0. 




1 



t °' 

o, 

V -1, 0, 4 
with Ti = 48 and for (mod 7 n ) the period is r 7 n = 48 ■ 7 n_1 



(52) 



16. To construct the matrix generator with eigenvalues in high fields GF[j¥p\ 
it is easier to use primitive polynomial of degree d over GF[^/p] the root of which 
coincides with the primitive element GF[Wp\. The primitive polynomial has the 



form El E2, 123 



d— 2 







(53) 



with coefficients /3i , /3 2 , • A* over GF[p]. The only problem is that this polynomial 
does not correspond to a matrix with unit determinant (2). But the last term {3d 
always can be represented as a power of the primitive element g of GF[p] (3 d = g k , 
therefore if we multiply the primitive polynomial (53) by A + g~ k we will get the 
polynomial which corresponds to a matrix with unit determinant 



(A + g- k )(X d + PiX + (3 2 \ d - 2 + ... + /3, 



X d+1 + (A + g- k )X d + ((3 2 + (3 1 g~ k )X d - 1 + ... + 1 = 0. 
To this polynomial corresponds the matrix generator of the form (33) 



(54) 



/ 0, 1, 0, 
0, 0, 1, 



.4 



d+l 



o, 
V -i, 



with period 



and on Z p n 



r p (A 



, \ 

., o 



, 1 

-({3 2 + {3 ig - k ), -((3 1 +g- k ) J 



p d - 1 



n— 1 



r pn (A d+1 ) = (p d -l)p 
This is our main result with the largest period of order p d+n . 



(mod p). 



(55) 



(56) 



(57) 



17. The example of the primitive polynomial on GF[7] with d = 10 is A 10 + 
A 9 + A 8 + 3 = and (54) has the form A 11 - A 10 - A 9 - 2A 8 - 4A + 1 = therefore 
the matrix is 



/ 0, 
0. 



A 



ii 



1, 
0, 1 



o \ 





0,0,0, ...,0,0,0, 1 



V 



(mod 7) 



(58) 



-1,4,0, ...,0,2, 1,1 J 

with period 7 10 - 1 and (7 10 - l^" 1 on Z pn . 

It is also useful to have the list of primitive polynomials on GF[2] ||22|| . Tables 
with larger ranges of d are available for GF[2]. In particular contain tables for 



d < 101, in |26j for d < 169 and in [27|] for d < 1001 with the corresponding period 
2iooo _ i 



d 

d~- 
d~- 
d = 



24 
= 81 
= 97 
= 127 
159 
d= 165 
d= 167 



^24 | \23 i \22 



+ X M + X zz + A 17 + 1 = 



A 81 + A 4 + 1 : 
A 97 + A 6 + 1 = 
A 127 + A + 1 






= 



d 



A 159 + A 31 + 1 = 
A 165 + A 31 + A 30 + A + 1 = 



A 167 + A + 1 = 0. 



In this case one can directly construct the matrices of the form (33) because the free 
term is equal to one. For the last polynomial we have 



/0, 
0. 



A 



167 







o \ 





(mod 2) 



(59) 



0, 0,0,0,0,0,1 

V i, 1,0,0,0,0,1 ) 

with period t 2 (A 2 q) = 2 167 — 1. Direct check of the eigenvalues, shows that eigenval- 
ues are not on a unit circle, therefore the K conditions (2,3) are satisfied. We have 
checked that for all primitive polynomials on hand this conditions are satisfied, so 
one can use any of them. 
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